Non-monotoic fluctuation-induced interactions between dielectric slabs carrying 

charge disorder 



O 

(N 

< 
(N 



o 

c3 



T3 

o 
o 



> 

00 

m 
p 

od 
O 

o 



X 



Jalal Sarabadani, 1,2 Ali Naji, 3 David S. Dean, 4 Ron R. Horgan, 3 and Rudolf Podgornik 2 ' 5 

1 Department of Physics, University of Isfahan, Isfahan 81746, Iran 
2 Department of Theoretical Physics, J. Stefan Institute, SI- 1000 Ljubljana, Slovenia 
8 Department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences, 
University of Cambridge, Cambridge CBS OWA, United Kingdom 
1 Laboratoire de Physique Theorique (IRSAMC),Universite de Toulouse, UPS and CNRS, F-31062 Toulouse, France 
''Institute of Biophysics, School of Medicine and Department of Physics, 
Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia 

We investigate the effect of monopolar charge disorder on the classical fluctuation-induced inter- 
actions between randomly charged net-neutral dielectric slabs and discuss various generalizations 
of recent results (A. Naji et al, Phys. Rev. Lett. 104, 06060f (20f0)) to highly inhomogeneous 
dielectric systems with and without statistical disorder correlations. We shall focus on the specific 
case of two generally dissimilar plane-parallel slabs, which interact across vacuum or an arbitrary 
intervening dielectric medium. Monopolar charge disorder is considered to be present on the bound- 
ing surfaces and/or in the bulk of the slabs, may be in general quenched or annealed and may 
possess a finite lateral correlation length reflecting possible 'patchiness' of the random charge dis- 
tribution. In the case of quenched disorder, the bulk disorder is shown to give rise to an additive 
long-range contribution to the total force, which decays as the inverse distance between the slabs 
and may be attractive or repulsive depending on the dielectric constants of the slabs. By contrast, 
the force induced by annealed disorder in general combines with the underlying van der Waals 
forces in a non-additive fashion and the net force decays as an inverse cube law at large separations. 
We show however that in the case of two dissimilar slabs the net effect due to the interplay be- 
tween the disorder-induced and the pure van der Waals interactions can lead to a variety of unusual 
non-monotonic interaction profiles between the dielectric slabs. In particular, when the intervening 
medium has a larger dielectric constant than the two slabs, we find that the net interaction can 
become repulsive and exhibit a potential barrier, while the underlying van der Waals force is attrac- 
tive. On the contrary, when the intervening medium has a dielectric constant in between that of 
the two slabs, the net interaction can become attractive and exhibit a free energy minimum, while 
the pure van der Waals force is repulsive. Therefore, the charge disorder, if present, can drastically 
alter the effective interaction between net-neutral objects. 



I. INTRODUCTION 

One of the most persistently made assumptions in the 
standard theory of (bio)colloid stability is that of a uni- 
form charge distribution on macromolecular surfaces [l[ . 
There are nevertheless many instances where charge pat- 
terns on macromolecular surfaces are inhomogeneous, ex- 
hibiting even a highly disordered spatial distribution. 
Among most notable examples are surfactant-coated 
surfaces and random polyelectrolytes and polyam- 
pholytes [3j . The macromolecular charge pattern in these 
systems is often frozen, or quenched, meaning that the 
charge distribution does not evolve after the assembly or 
fabrication of the material. This is the case if the inter- 
action between the surface charge carriers and the sur- 
faces is of a chemical nature involving interactions that 
far exceed the thermal energy scale k^T. Alternatively 
charge distributions can exhibit various degrees of an- 
nealing when interacting with other macromolecules in 
aqueous solutions as is, for example, the case for charge 
regulation of contact surfaces bearing weak acidic groups 
in aqueous solutions 

Recently it has also been realized that ultrahigh sen- 
sitivity experiments on Casimir (zero temperature and 
ideally polarizable surfaces) and van der Waals (finite 



temperature and non-ideally polarizable surfaces) inter- 
actions between surfaces in vacuo 0, Q can be prop- 
erly understood only if one takes into account the dis- 
ordered nature of charges on and within the interacting 
surfaces Possible causes of sample- and history- 

specific charge disorder in this case include the patch ef- 
fect, where the variation of the local crystallographic axes 
of the exposed surface of a clean polycrystalline sample 
can lead to a variation of the local surface potential Q . 
Amorphous films deposited on crystalline substrates can 
also show a similar type of surface charge disorder, show- 
ing a grain structure of dimensions sometimes larger than 
the thickness of the deposited surface film [l(|. On the 
other hand, adsorption of various contaminants can influ- 
ence the nature and type of the surface charge disorder. 

Since the nature and distribution of the charge dis- 
order in any of the force experiments is in general sel- 
dom known, we pursue a strategy of assessing the con- 
sequences of different a priori models of the distribution 
of charge disorder. In what follows we will thus assume 
that the charge disorder stems from randomly distributed 
monopolar charges which may be present both in the bulk 
and/or on the interacting surfaces and can be either an- 
nealed or quenched. In the quenched case, the disor- 
der charges are frozen, whereas in the annealed case, the 
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disorder charges are subject to thermal fluctuations at 
ambient temperature and can thus adapt themselves in 
order to minimize the free energy of the system. We 
do not deal with effects due to disorder in the dielectric 
response of the interacting media [111 ], which presents 
an additional source of disorder meriting further study. 
We have shown in our previous works [H, [l2l [l3j that 
the type and the nature of the charge disorder induces 
marked changes in the properties of the total interaction 
between apposed bodies and the force arising from the 
presence of disordered charges can dominate the under- 
lying Casimir-van der Waals (vdW) effect at sufficiently 
large separations. The analysis of the interaction finger- 
print of the charge disorder can be useful in assessing 
whether the experimentally observed interactions can be 
interpreted in terms of disorder effects or are due to pure 
Casimir-vdW interactions. In this study we will continue 
with the assessment of this interaction fingerprint in the 
case two semi-infinite net-neutral dielectric slabs (sepa- 
rated by a layer of vacuum or an arbitrary unstructured 
dielectric material as shown in Fig. [I}, by generalizing 
our formalism to include two different mechanisms: 

i) the "patchy" distribution of the charge disorder, 
which may arise due to finite correlations between 
the carriers of the charge disorder. This is a rather 
realistic assumption in view of the graininess ob- 
served in interacting surfaces [1, [T3| • 

ii) the dielectric inhomogeneity effects in highly asym- 
metric systems, where the two slabs and the inter- 
vening medium may in general have different di- 
electric properties. 

For interacting bodies carrying no net charge, the 
part of the interaction due to monopolar quenched 
charge disorder can be seen from a full quantum field- 
theoretical formalism [l4| to be coupled directly to the 
zero-frequency or classical vdW interaction (correspond- 
ing to the zero-frequency Matsubara modes of the elec- 
tromagnetic field). The situation would be more compli- 
cated in the case of annealed disorder if one is to account 
for the higher-order Matsubara frequencies and a general 
treatment is missing at present. In this work, we shall 
examine the effects of quenched and annealed charge dis- 
order by focusing specifically on the zero-frequency vdW 
interaction. It is well known that the higher frequency 
Matsubara terms in the total vdW interaction free en- 
ergy become relatively unimportant at sufficiently high 
temperatures and/or sufficiently large intersurface sep- 
arations In this paper we will restrict ourselves to 
this regime which is also most relevant experimentally 
Q and consider the short-distance quantum effects else- 
where 

As we shall demonstrate, the interplay and competi- 
tion of the underlying vdW effect between the two semi- 
infinite slabs and the interaction induced by the charge 
disorder can give rise to a variety of different interaction 
profiles, including, most notably, non-monotonic interac- 



tions as a function of the inter-slab distance. In partic- 
ular, when the two dielectrics interact across a medium 
of higher dielectric constant, the disorder force can be- 
come strongly repulsive and thus generate a potential 
barrier when combined with the vdW force, which is at- 
tractive in that case. While, in the situation where the 
intervening medium has a dielectric constant in between 
that of the two slabs [Hj], the disorder can generate an 
attractive force, which can balance the repulsive vdW 
force, leading thus to a stable bound state between the 
two dielectric slabs. Therefore, the charge disorder, if 
present, can provide an intrinsic mechanism to stabilize 
interactions between net-neutral bodies. These features 
emerge because the disorder-induced interactions typi- 
cally decay more weakly with the separation than the 
vdW interaction and depend strongly on the dielectric 
inhomogeneities in the system. 

The above features vary depending on whether the dis- 
order is assumed to be quenched or annealed. The ef- 
fects due to annealed charge disorder can be distinctly 
different from those of the quenched charge disorder. In 
the case of annealed charges, the effective interactions 
turn out to decay more rapidly and in fact in a similar 
way as the pure vdW interaction, whereas in the case of 
quenched charges the net interactions decay more weakly 
with the separation, i.e., as ~ D~ a where a = 1 if the 
disorder is present in the bulk of the slabs and a = 2 if 
the disorder is confined to the bounding surfaces of the 
slabs. 

The organization of the paper is as follows: In Section 
HI] we introduce the details of our model and the formal- 
ism used in our study. The general results are derived in 
Sections IIII Al and IIII Bl for the case of quenched and an- 
nealed disorder, respectively. We shall proceed with the 
analysis of the effective interaction between two identical 
slabs in Section IIV1 where we shall focus mainly on the 
effects due to spatial correlation in the distribution of 
the charge disorder. In Section [Vj we study the interac- 
tion in dielectrically asymmetric systems by considering 
two dissimilar dielectric slabs that interact across an ar- 
bitrary dielectric medium. The results and limitations of 
our study are summarized in Section IVII 



II. MODEL 

Let us consider two semi-infinite slabs of dielectric con- 
stants £1 and £2 and temperature T with parallel planar 
inner surfaces (of infinite area S) located normal to the 
z axis at z = ±D/2 (see Fig. [T]). The inner gap is filled 
with a material of dielectric constant e m . The inhomo- 
geneous dielectric constant profile for this system is thus 
given by 

f ex z< -D/2, 
e(r) = { e rn \z\ < D/2, (1) 
[ e 2 z > D/2. 
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FIG. 1: (Color online) We consider two semi-infinite net- 
neutral slabs (half-spaces) of dielectric constant si and £2 
interacting across a medium of dielectric constant e m . The 
monopolar charge disorder (shown schematically by small 
light and dark patches) is distributed as random patches of 
finite typical size (correlation length) in a layered structure 
in the bulk of the slabs and on the two bounding surfaces at 
z = ±D/2. It may be either quenched or annealed. 



We shall assume that the two dielectric slabs have a dis- 
ordered monopolar charge distribution, p(r), which may 
arise from randomly distributed charges residing on the 
bounding surfaces [p s (r)] and/or in the bulk [pb(r)], i.e., 
p(r) = p s {r) + p b (r). The charge disorder is assumed to 
be distributed according to a Gaussian weight with zero 
mean (i.e., the slabs are net neutral), and the two-point 
correlation function 

({p(r)p(v')))=g(e-Q';z)S(z-z'), (2) 



and that the surface and bulk variances are given by 

g s (z) = e 2 [g ls S(z + D/2) + g 2s 8{z - D/2)], (6) 

( gi b e% z < -D/2, 
9b(z) = < \z\ < D/2, (7) 

[ g 2b e\ z > D/2. 

The lateral correlation between two given points is typ- 
ically expected to decay with their separation over a finite 
correlation length ("patch size"), which could in general 
be highly material or sample specific. However, the main 
aspects of the patchy structure of the disorder can be in- 
vestigated by assuming simple generic models with, for 
instance, a Gaussian or an exponentially decaying cor- 
relation function. Without loss of generality, we shall 
choose an exponentially decaying correlation function ac- 
cording to the two-dimensional Yukawa form 

c iQ (x) = -l r K fM N ) ) ( 8 ) 

where £i a represents the correlation length for the bulk or 
surface disorder (a = b, s) in the i-th slab (i = 1,2). The 
case of a completely uncorrelated disorder [8] follows as a 
special case for £ ia — > from our formalism. We should 
emphasize that the correlation assumed for the bulk dis- 
order is only present in the plane of the slab surfaces and 
not in the direction perpendicular to them. This assump- 
tion is wholly justified only for layered materials. In all 
the other cases of the bulk disorder one would normally 
expect the same correlation length in the direction per- 
pendicular to the surfaces. We will deal with this model 
in a separate publication. 



where ((■••)) denotes the average over all realizations of 
the charge disorder distribution, p(r). We have thus as- 
sumed that there are no spatial correlations in the per- 
pendicular direction, z, while, in the lateral directions 
q = (x,y) (in the plane of the dielectrics), we have a fi- 
nite statistically invariant correlation function whose spe- 
cific form may depend on z as well. This implies that 
the charge disorder is distributed in general as random 
"patches" in a layered structure in the bulk of the slabs 
as well as on the bounding surfaces. 

The total correlation function can be written as the 
sum of the surface (s) and bulk (b) contributions 

G{q-q';z) = g s {z)c s (g- Q';z)+g b (z)c b {g- q';z). (3) 

For the slab geometry, we generally assume that the lat- 
eral correlation functions may be different for the two 
slabs, i.e. 



c s (x;z) 
Cb(x; z) 



Cl s (x) 
C2s(x) 



-D/2. 
D/2, 



ci b (x) z < -D/2, 
\z\<D/2, 
c 2b (x) z > D/2, 



(4) 
(5) 



III. FORMALISM 

The partition function for the classical vdW interac- 
tion (the zero-frequency Matsubara modes of the elec- 
tromagnetic field) may be written as a functional integral 
over the scalar field </>(r), 



Z[p{r)} = / [Xty(r)] e-^^'lfWl, 



(9) 



with f3 = l/k-eT and the effective action 
5[0(r);p(r)]= f dr [| eo e(r) (V0(r)) 2 + i p(r)0(r)] . 



(10) 

The above partition function can be used to evaluate 
averaged quantities such as the effective interaction be- 
tween the dielectric bodies. However, since the charge 
distribution, p(r), is disordered, it is necessary to aver- 
age the partition function over different realizations of 
the charge distribution. The averaging procedure differs 
depending on the nature of the disorder. In what fol- 
lows, we consider two idealized cases of either completely 
quenched or completely annealed disorder fl6j (the inter- 
mediate case of partially annealed disorder is also ana- 
lytically tractable [l3[ but will not be considered here). 
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A. Correlated quenched disorder 

The quenched disorder corresponds to the situation 
where the disorder charges are frozen and can not fluc- 
tuate and equilibrate with other degrees of freedom in 
the system. As it is well known, the disorder average 
in this case must be taken over the sample free energy, 
lnZ[p(r)], in order to evaluate the averaged quantities 
[lt| . Therefore, the free energy of the quenched model is 
given by 



quenched 



= -((lnZ[p(r) 



(11) 



The Gaussian integral in Eq. © as well as the disorder 
average can be evaluated straightforwardly in this case 
yielding 



Quenched = iTrlnG'- 1 (r,r') 



(12) 



J dr dr' Q{q- g'; z)5(z - z')G(r, r'), 



where G(r, r') is the Green's function defined via 

eoV • [e(r)VG(r, r')] = -S(v - r'). (13) 

Note that the free energy (|12p of the quenched model is 
expressed in terms of different additive contributions, i.e. 



•J quenched — J~ vdW ~t~ ^~b J~s 



(14) 



The first term above is nothing but the usual contribution 
from the zero-frequency vdW interaction, 



/37"vdW = i r tVlnG- 1 (r,r') 



(15) 



which is always present between neutral dielectrics even 
in the absence of any monopolar charge disorder. The 
second and the third terms represent contributions from 
the bulk and surface disorder (a = b, s) 

= | J "dr dr' g a (z)c a (g- g';z)S(z - z')G(v,v'). 

.(16) 

The quenched expression (|12p is valid for any arbitrary 
disorder correlation function Q(g — g';z) and dielectric 
constant profile e(r). In what follows, we shall focus on 
the particular case of planar dielectrics by employing Eqs. 
©, ©-©• In this case, the zero-frequency vdW contri- 
bution is obtained (per k B T and unit area) as Q 



vdW 



IIW M1 - A ^ W) - (17) 



The force associated with this contribution, / v( jw 
—dFvdw/dD, follows as 



Pfvdw Li 3 (A!A 2 ) 



S 8irD 3 ■ 

The dielectric jump parameters are defined as 

a £i £m 



Si + £r, 



(18) 



(19) 



at each of the bounding surfaces (i = 1,2), and Li3( 
the trilogarithm function defined by 



\ -> z 



(20) 



The dielectric constants of the system enter the expres- 
sion p7|) via the terms Ax and A 2 . If AiA 2 > then 
the vdW interaction between the slabs is attractive, oth- 
erwise the vdW force is repulsive. 

The expression in Eq. (TTIj)) can be calculated in a 
straightforward manner. We obtain the bulk disorder 
contribution as 



S 

dQ 



-£ B £ m x 

P -2QD 



(21) 



Q l-AiAze" 2 ^ 



9lbC\b{Q) * , 92bC2b{Q) 



(ei+£r. 



:A 



Ai 



and the surface disorder contribution as 



(3F S 
S 



= -2£ B £ m x 

P -2QD 



(22) 



1- AiA 2 e- 2C ? D 
at all separations D, where 

£ B = (3e 2 /(<iTr£ 



guCi s (Q) ^ | g 2s c 2s (Q) 
(ei+e m ) 2 2 (e 2 +e m ) 2 



Ai 



(23) 



is the Bjerrum length in vacuum (£ B ~ 56.8 nm at room 
temperature), and Ci a (Q) is the Fourier transform of the 
correlation function Cj Q (x). For the particular form of 
Cj Q (x) chosen in Eq. ([5)1, we have a Lorentzian Fourier 
transform (for £ = 1,2 and a = b, s) as 



Cia (Q) 



CQ 2 + i 



(24) 



The total force between the dielectric slabs carrying 
quenched charge disorder thus follows from / qU enchcd = 
— <9J-q U cnchcd/<9-D and the preceding free energy expres- 
sions as 



/3/< 



quenched 

dQ e 



Li 3 (A!A 2 ) 



-2QD 



87TD 3 

gucib(Q) 



2£ B £r, 



(25) 



o (l-AiAae" 2 ^ 



(£i+£ m y 



2Qgi s cis(Q) 
(£l+£ m ) 2 



A 2 



(e 2 + e. , 
2Q.g 2s c 2s (Q) 

(£2+£m) 2 



A! 



We notice that the quenched disorder contribution to 
the force has a number of interesting features. We see 
that the forces generated due to the charge disorder in 
slab 1 and slab 2 are independent and additive. If we 
remove the charge disorder from either slab we see there 
is still a contribution coming from the other slab. This is 
because the charge distributions on average do not inter- 
act with each other. As the charge distribution in oppos- 
ing slabs is uncorrelated, the average force on a charge 
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in slab 1 due to a charge in slab 2 is zero as the charge 
in slab two is equally likely to have a positive charge 
as a negative charge. In fact the charges in slab 1 are 
only correlated with their image charges in slab 2, this 
means that on average the charge distribution in slab 1 
only interacts with its image in slab 2 and vice versa. 
The sign of the interaction between charges in slab 1 and 
their images in slab 2 depends on A2. If A2 is positive, 
that is the dielectric constant of slab 2 is greater than 
that of the intervening material, £2 > e m , then the force 
due to the charge in slab 1 is attractive. Therefore, the 
contribution of the charge in each slab to the net force 
depends on the dielectric contrast between the opposing 
slab and the intervening medium. Thus, while the vdW 
contribution depends on the product A1A2, the contri- 
bution from the quenched charge distribution in slab 1 
depends on variance <7i s ,& of the charge on the surface or 
bulk of this slab but on A2 of the slab 2. This leads to 
a rich phenomenology of possible net interactions. For 
example, if Ai and A 2 are both positive, the vdW force 
is attractive and so are the forces due to both charge 
distributions. However, if Ai and A2 are negative the 
vdW force is again attractive but the force due to both 
quenched charge distributions is repulsive. Also, note 
that while the vdW term shows a standard 1/D 3 decay 
with the surface separation D, the disorder contribution 
turns out to have a much weaker decay with the sepa- 
ration Q as we shall analyze further in the forthcoming 
sections. Therefore, the interplay between these differ- 
ent contributions can lead to characteristically different 
interaction profiles depending on the system parameters. 

B. Correlated annealed disorder 

In the annealed model, disorder charges are assumed 
to fluctuate in thermal equilibrium with the rest of the 
system, and thus in particular the charge distribution in 
the two slabs can adapt itself to minimize the free en- 
ergy of the system. In this case, the disorder degrees of 
freedom should be treated statistically on the same foot- 
ing as other degrees of freedom, which implies that the 
disorder average must be taken over the sample parti- 
tion function Z[p(r)] [l6j]. Hence, the free energy of the 
annealed model reads 

^annealed-- In ((Z[p(r)])). (26) 

The annealed free energy may be evaluated using Eqs. 
©, © and (HU) as 

Annealed - ^Tr In p- 1 (r, v')+pg(Q-Q'- z)S(z - z')) . 

(27) 

Note that unlike the quenched case, Eq. (fl"2"j) . the disor- 
der contributions cannot be in general separated from the 
pure vdW contribution when the disorder is annealed. 

In the case of two interacting planar dielectrics, we 
shall make use of Eqs. (H} and ©-(0 in order to cal- 
culate the fluctuational trace-log of the modified inverse 



Green's function G 1 (r, r') + j3Q(g— g'; z)S(z-z') in Eq. 
(|27p . By employing standard procedures [17] , we find 

fal = I/ i^ln[l-A lg (g)A 29 (Q) e - 2 ^], 

(28) 

where we have (for i = 1,2) 

A »fl(<3) = (29) 

_ £ m Q-Si^/Q 2 +^t B g. lb c. lb (Q)/~i-Ant B g ls c is (Q) 
e m Q+Ei v /Q 2 +47r^B5ib c ib (Q)/si + ^t^gis <ks (Q) 

The total force, /annealed = -dJanncaicd/d-D, in the an- 
nealed model thus follows as 

/3/annealed 1 [p^n A lg (Q) A 2g (Q) e' 2 ^ 

S 2nJ * ^ 1 - A lg (Q)A 2g (Q) e~^ D ' 

The annealed interaction free energy, or the correspond- 
ing force, has a form reminiscent of the zero- frequency 
vdW interaction between two semi-infinite slabs [l8j with 
renormalized dielectric missmatch, Aj S (Q), depending 
on the in-plane wave-vector Q. Indeed, one can ar- 
gue for the following exact correspondence based on 
the vdW interactions between media with volume and 
surface embedded mobile charges [TH, Ixsj : the term 

^/Q 2 +47r^B.9ifc Cib(Q)/ei corresponds to bulk Debye-like 
"screening" , stemming from the annealed bulk disorder 
response to local electrostatic fields. It can be obtained 
alternatively by analyzing the fluctuational modes of the 
Debye-Hiickel fields as opposed to the Laplace fields [l8[ . 
The term 47r^B<?is Ci S {Q), on the other hand, stemms from 
the response of the mobile charges confined to the di- 
electric surfaces to electrostatic fields, corresponding to 
surface disorder-generated "screening" [lj|. It can be 
derived alternatively by considering the response of the 
surface polarization to local electrostatic fileds 19]. The 
combination of the two terms in Eq. (1301) can not be 
written as a sum of two terms depending linearly on Ci b 
and Ci S . This signifies that the bulk and surface disorder 
effects for the annealed case are obviously not additive 
and can not be analyzed separately. 

IV. SYMMETRIC CASE OF TWO IDENTICAL 
SLABS 

We now analyze the preceding analytical results in the 
situation where the two dielectric slabs are identical, i.e., 
we have e, = e p (Aj = A), g ls = g 8 , g ib = g b , and 
Cis = db = 6 for both slabs i = 1,2. 

A. Quenched disorder-induced interactions 

In the case of quenched disorder, one can expand the 
Lorentzian correlation function c ia (Q), Eq. (|24p . in pow- 
ers of the dimensionless ratio of the correlation length to 
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FIG. 2: a) Magnitude of the rescaled total force, f3\f\l B /S (Eq. I|25|)). between two identical net-neutral dielectric slabs in 
vacuum (e m = 1) bearing quenched monopolar charge disorder as a function of the rescaled distance, D/Ib- The results are 
plotted here for fixed e v = 10, g s — (no surface disorder) and bulk disorder variance gi, = 5 x 10~ 8 nm~ 3 and varying disorder 
correlation length (,/Ib = 0, 200, 10 3 , 10 4 , 10 s (from top to bottom). Inset shows the ratio of the total force to the pure zero- 
frequency vdW force f|34f) between the slabs in the absence of charge disorder for the same range of D. b) Same as (a) but here 
we fix e p = 10, g s = 0, £/Ib = 200 and vary the bulk disorder variance in the range gb = 10 -7 , 10 -8 , 10 -9 , 10 -10 , 10~ n nm~ 3 
(from top to bottom). Solid curve is the pure vdW force (|34|) . c) Same as (b) but here we fix £/Zb = 200, gt = 5 x 10~ 8 nm -3 , 
g s = g^ 3 and vary the dielectric constant of slabs as e p = 2 (top dotted curve), 10, 40, 100 (bottom dot-dashed curve). All 
graphs are plotted in log-log scale. 



the intersurface separation, £/D, and thus express the 
total force (|25|) in the form of a series expansion as 



/quenched f 



(n) 



(31) 



where we have 
/3/ (0) _ qUb A 2£ m5 ^ B ln(|l-A 2 | 



Li 3 (A 2 ) 



S 



2e p D 



( £m + £p ) 2 AD 2 



8ttD 3 ' 

corresponding to the free energy of the system in the 
presence of a completely uncorrelated disorder (£ = 0), 
and 



/3/ ( 



S (e 



°° tin 



r ' n—1 



C'2n+lLi2n(A 



2/i 



(-l) n ^C2„+ 2 Li 2 „ + i(A J ), (33) 



(e m + £ P ) 2 A£2 ^ 

v ^ 7 n— 1 

with C„ = r(n)/2™. The above series expansion is most 
suitable for the situation where £/D is small. 

Let us first consider the leading-order uncorrelated dis- 
order free energy /(°) in Eq. (|32|) . This contribution 
exhibits a few remarkable features [8]. First, it shows 
a sequence of scaling behaviors with the separation that 
stem from different origins: a leading 1/D term due to 
the quenched bulk disorder, a subleading 1/D 2 term from 
the surface charge disorder, and the pure vdW term that 
goes as 1/D 3 , i.e. 



PfvdW 



Li 3 (A 2 ) 
8ttD 3 ' 



(34) 



While the vdW term is always attractive in a symmet- 
ric system, the disorder contributions (first and second 
terms in Eq. (|32p ) are attractive when the dielectric mis- 
match A > (i.e., when the dielectric constant of the 
intervening medium is smaller than that of the slabs) 
and repulsive otherwise. The former case (A > 0) has 
been investigated in detail previously in the context of 
two identical dielectric slabs in vacuum (e m = 1) 



Note that since the forces induced by an uncorrelated 
disorder exhibit a much weaker decay with the separa- 
tion, they may completely mask the standard vdW force 
(which is always present regardless of any charge dis- 
order) at sufficiently large separations One might 
expect that globally electroneutral slabs would exhibit a 
dipolar-like interaction force to the leading order rather 
than the monopolar forms 1/D (or 1/D 2 ) obtained for 
the bulk (or surface) charge distribution. As noted be- 
fore, the physics involved is indeed subtle as the disorder 
terms result from the self-interaction of the charges with 
their images (which follows from G(r,r), Eq. (fl"2"j) . in 
the limit of zero correlation length) and not from dipolar 
interactions (which come from an expansion of G(r,r') 
when |r — r'| is large). Statistically speaking each charge 
on average (as any other charge has an equal probability 
of being of the same or opposite sign) only sees its im- 
age, thus explaining the leading monopolar form in the 
net force. 



The next leading correction due to disorder correla- 
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tions for small £/D follows from Eq. (|33[) as 

S ~ (e p + e m ) 2 D 3 + (e m +e p ) 2 AD* ' 

(35) 

which has an opposite sign to the zeroth-order term, and 
thus tends to weaken the zeroth-order effect. Note how- 
ever that the higher-order corrections in Eq. (|33|) alter- 
natively change sign and their net effect leads to a total 
force which can be evaluated numerically, e.g., directly 
from Eq. (1251) . The results are shown in Fig. [2] for the 
case of two identical dielectric slabs in vacuum. As seen 
in Fig. EK, the attractive disorder-indued forces indeed 
diminish (albeit rather slowly as seen in the inset) as the 
disorder correlation length is increased (by several orders 
of magnitude from £ = up to £ ~ 5mm in actual units 
in the figure). For a large correlation length or at small 
separation (large £/D), the total force thus tends to the 
non-disoredred vdW force ([34]) that scales as 1 /D 3 (thick 
solid line) , while for a small correlation length or at large 
separation (small £/D), the force increases and shows 
a crossover to the maximal uncorrelated disorder value 
([32]) that scales as 1/D (top dotted line). The latter is 
obviously from the bulk disorder which, if present, gives 
rise to the most dominant effects. 

The crossover behavior depends also significantly on 
the disorder variance and the dielectric mismatch (Figs. 
[2]b and c). Similar trends as described above can be 
observed by varying the bulk disorder variance, gb, as 
shown in Fig. [2Jd (here gb varies within the typical range 
1CP 11 — 10~ 7 nm~ 3 corresponding to impurity charge den- 
sities of 10 10 - 10 14 e /cm 3 [2(ij). The results in Fig. Hfc 
however reveal a nonmonotonic dependence on the di- 
electric mismatch: The disorder contribution tends to 
decrease both for small and large slab dielectric constant. 
This can be seen directly from Eqs. (|32|) and ([33]) as the 
disorder-induced force vanishes in the limits e p — > oo or 
A — » 1 (perfect conductor) and A — > (homogeneous 
medium, in which case the vdW force vanishes as well). 

It is interesting to note that the first-order corrections 
due to bulk correlations scales as which is similar to 

the vdW contribution. Hence, these corrections tend to 
renormalize the ideal zero- frequency Hamaker coefficient, 
A = 3(fc B T)Li 3 (A 2 )/4, associated with the vdW term 
(defined via f v dw/S = —A/ (6nD 3 )) to an effective value 
given by 



A cS = ^Li 3 (A 2 ) 



3ff fe e 2 e m A£ 2 Li 2 (A 2 
2e (e p + e m ) 2 
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FIG. 3: a) Ratio of the total force (|30|l to the zero-frequency 
vdW force (|34[1 between two identical net-neutral dielectric 
slabs in vacuum (s m = 1) bearing annealed monopolar charge 
disorder as a function of the rescaled distance, D/Ib- The 
results are plotted here for fixed e p = 10, g s — (no sur- 
face disorder), gt, = 5 x 10~ 8 nm -3 and varying disorder 
correlation length f/i B = 0, 200, 10 3 , 10 4 (from top to bot- 
tom). Annealed curves are bounded by the perfect conduc- 
tor result, Eq. (|37[) (top solid line), for large disorder and 
the vdW result, Eq. (|34[) (bottom solid line), for no dis- 
order, b) Same as (a) but here we have fixed e v = 10, 
g 3 — 0, £/Zb = 200 and bulk disorder variance varied in 
the range g b = 10" 7 , 10" 8 , 10 -9 , 10 -10 , 10~ n nm~ 3 (from 
top to bottom), c) Same as (a) but here we fix g 3 = 0, 
gt, = 5 X 10~ 8 nm 
of the slabs as e ; 
closer view of the curves for e p = 40, 100 (from top). Top 
dotted lines correspond to Eq. (|37[) . All graphs are plotted 
in log-log scale. 



-10 

Same as (a) but 
~ 3 ,£ = 0, and vary the dielectric constant 
p — 2,10,40,100 (from top). Inset shows a 



This value is larger (smaller) than the ideal value A when 
A > (A < 0) and can even change sign. 



B. Annealed disorder-induced interactions 

In the case of annealed disorder, the total force between 
the dielectric slabs is given by Eq. (|3T)]) . It shows that the 



force in this case cannot be expressed simply as a sum of 
different additive contributions as in the quenched case. 
Rather, it follows from the combined effect of the an- 
nealed disorder fluctuations and the underlying vdW ef- 
fect, which again reflects the fact that the disorder charge 
distribution in this case can adapt itself in order to min- 
imize the total free energy of the system. 
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FIG. 5: The rescaled magnitude of the (attractive) force, 
Eq. (|30[) . between two identical net-neutral dielectric slabs 
bearing annealed monopolar charge disorder in a dielectri- 
cally homogeneous system (ei = £2 = Em) as a function 
of the rescaled distance, D/Ib- Here we plot the resulting 
force for ei = £2 — £m — 1, 10, 40, 100 (from top to bottom) 
and for uncorrelated bulk disorder in both slabs with £ = 0, 
gb = 5 x 10 -8 nm~ 3 and g s — 0. Top solid line shows the 
universal limiting expression ()37|) . 



FIG. 4: a) Same as Fig. [3^ but here we fix e v = 50, g s = 0, 
<?b = 5 x 10 -8 nm -3 , £ = 0, and vary the dielectric constant of 
the intervening medium as e m = 2, 10, 20, 40 (from bottom to 
top), b) Same as (a) but here we fix the dielectric constant of 
the slab as e p = 10 and vary that of the intervening medium 
as Em = 30, 40, 60, 80 (from top to bottom). Top dotted lines 
correspond to Eq. ((37)) . 



The resulting total force is shown in Fig. [3] for the 
symmetric case of two identical dielectric slabs in vac- 
uum. For the sake of presentation, we have plotted the 
ratio of the annealed force (150]) to the vdW force (IMl) . 
which clearly shows that the net force is bounded by two 
limiting laws, namely, the ideal limiting force 

/ideal _ C(3) ,„_,. 
S ~ 8tt£> 3 ' [6/) 

which is similar to the one obtained between perfect con- 
ductors @, and the vdW force 433} 

in the absence of 

any disorder charges. These constitute the upper and 
the lower bounds for the annealed force in the general 
symmetric case with A > 0, i.e. 

|/vdw| < I /annealed I < | /ideal |- (38) 

Note that the above limiting values can be established 
systematically from the general expression and are 
strictly attractive (negative sign). The lower bound vdW 
force is non-universal (i.e., material dependent) and is 
obtained asymptotically in the limit of small disorder 
variance (gb, g s 0) or small separationss. The upper 
bound is universal and follows in the limit of strong dis- 
order variance (gb or g s — > 00) or large separations. The 
crossover from one limit to the other can be achieved by 
increasing the distance or the disorder variance as seen 
in Fig. El 



In the annealed case, the disorder correlations play a 
similar role as in the quenched case and tend to weaken 
the disorder-induced forces. As seen in Fig. [3k,, the mag- 
nitude of the force drops from its value for an uncorre- 
lated disorder (£ = 0), and also the crossover between 
the two limiting behaviors as discussed above is 'delayed' 
when the disorder has a finite correlation length, £. 

It is also remarkable to note that Eq. (|3"T)) is obtained 
in general as the large-distance (D — > 00) behavior for 
the total force between any two arbitrary dielectric slabs 
regardless of their dielectric constant and disorder vari- 
ance. It thus demonstrates the intuitive fact that dielec- 
tric slabs with annealed charges tend to behave asymp- 
totically in a way similar to perfect conductors. This is 
one of the distinctive features of the annealed disorder 
as compared with the quenched disorder, whose effects 
depend significantly on the material properties. The de- 
viations due to material properties and the disorder vari- 
ance in the annealed case contribute a repulsive sublead- 
ing force in the symmetric case, which can be determined 
from a series expansion at large separations. If the sur- 
face disorder is negligible as compared with the bulk dis- 
order (g s <C gb^B), we obtain (up to the first few leading 
orders) 

/9/annealed _ C(3) 3g m C(3) _ 3 £ ^C(3) 

S ~ 8ttD 3 v /64 7 r 3 (? b £B£ P £' 4 8TT 2 g b e B epD 5 
15£ 2 e m C(5) 45^4^(5) m 

2 v / 6ATT 3 g b £BSpD e 16ir 2 g b eBe p D 7 ' 

While in the opposite situation where the bulk disorder 
is negligible (g s ^> gb£ B ), we obtain 

/9/annealed C(3) 3 £m C(3) 15g m g 2 C(5) 

S ~ 8irD 3 16TT 2 g s £ B D 4 l6ir 2 £ B g s D 6 ' 

(40) 
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We should emphasize that, although the annealed force 
decays at large separations in a similar fashion as the 
pure vdW force, its magnitude can nevertheless exceed 
the vdW force by a few orders of magnitude if the dielec- 
tric constant of the intervening material s m is increased 
toward that of the slabs (see Fig. |4ji). The asymptotic 
behavior of the annealed force is similarly described by 
the pure vdW result and the ideal expression ([57)1 
even if the dielectric constant of the (identical) slabs is 
smaller than that of the intervening medium, A < 0. 
However, in this case, these two limiting results do not 
constitute the upper and lower bound limits for the total 
force as the force ratio / / / v dW exhibits a nonmonotonic 
behavior with a local minimum at some intermediate sep- 
aration between the slabs (Fig. |4Jd). 

The special case of a dielectrically homogeneous sys- 
tem, £\ =62 = e m , constitutes another example where 
annealed and quenched disorder effects differ on a quali- 
tative level and may thus be easily distinguished. In this 
case, the total force ([25]) due to quenched disorder van- 
ishes trivially at all separations, while the total force due 
to annealed disorder remains finite. Since the vdW con- 
tribution vanishes in a dielectrically homogeneous system 
as well, the total force in this case comes purely from the 
electrostatic interactions of annealed charges in the two 
slabs (Fig. EJ). 



V. ASYMMETRIC CASE OF TWO DISSIMILAR 
SLABS 

So far we have considered only the case of a symmet- 
ric system composed of two identical semi-infinite slabs. 
In practice, however, one may often deal with a situa- 
tion where the dielectric constant or disorder variance 
of the two slabs are different. In this case, the result- 
ing fluctuation-induced interactions may exhibit qualita- 
tively different features as compared with the fully sym- 
metric case that we shall explore further in this Section. 



A. Interaction between a disordered and a 
disorder-free slab 

Let us first consider briefly the situation where the two 
slabs are dielectrically identical, E\ = £2 (Ai = A2 = 
A), but bear different degrees of quenched or annealed 
monopolar charge disorder. For the sake of simplicity, let 
us assume that one slab is disorder free (gib,gis = and 
£1 =0), whereas the other slab contains disorder charges 
of bulk and surface variance 32b and g2 S , and correlation 
lengths ^2fc = 62s = £2 (recall that in any case the net 
monopolar charge in each slab is taken to be zero). 

This case is particularly interesting because it shows 
the interconnection between the disorder and the dielec- 
tric inhomogeneity in the system. In the quenched case, 
we find that the disorder contribution to the total force 
is nonzero (even though one of the slabs does not con- 
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FIG. 6: Same as Fig. [3^l but for two dissimilar dielectric 
slabs in vacuum (e m = 1) with one slab being disorder free 
(gib, gis = 0, £1 = 0) and the other slab containing annealed 
charge disorder of variance g2b = 5 x 10 -8 nm -3 , g2s = 0. 
Here we fix e p = 10 and vary the disorder correlation length 
as &/Ib = 0, 200, 10 3 , 10 4 (from top to bottom). The results 
in this case are bounded by the limiting values given by Eqs. 
fl34[) and J42J) (solid lines labeled by Li 3 (A 2 ) and Lis (A), re- 
spectively); see Eq. (|4ip . The top solid line (labeled by £(3)) 
is from Eq. ([37]) . 



tain any disorder charges), and is simply given by half 
the value obtained in the fully symmetric case (see Eq. 
(|3ip ) in the previous section if we set g2b = 36, 32s = 3s 
and £2 = £ (the vdW contribution is the same in both 
cases). This result follows straightforwardly from the 
general quenched expressions (|25|) or (|Alj) and again re- 
flects the fact that the quenched contribution is basically 
due to the interaction of disorder charges with their im- 
age charges (as these are the only 'charges' with which 
they are 'correlated' across the intervening gap). This 
is why the disorder forces in the quenched case depend 
essentially on the dielectric jump across the bounding 
surfaces as discussed before. 

In the annealed case, the disordered dielectric slab 
tends to behave asymptotically as a perfect conductor, 
while the disorder-free slab behaves as a dielectric mate- 
rial. This lowers the magnitude of the maximum net force 
that can be achieved in this system as compared with the 
case where both slabs bear annealed charges; the latter 
case is obviously more favorable thermodynamically as 
the system can achieve a lower free energy. The net an- 
nealed force in vacuum still falls between two well-defined 
limits as shown in Fig. [51 The two bounding limits are 
given by 



|/vdw| < |/i 



annealed 



<l/* 



(41) 



where / v dW, Eq. (|34[) . is obtained in the limit of weak 
disorder or at small separations, and /*, defined as 



S 



Li 3 (A) 
' 8tt£> 3 ' 



(42) 



is obtained in the limit of strong disorder or large sep- 
aration. Thus, in marked contrast with the fully sym- 
metric case, we find that the large-distance interaction 
here becomes non-universal and could be either attrac- 
tive (A > 0) or repulsive (A < 0). This points to the 
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FIG. 7: a) The rescaled total force, PfV^/S (Eq. I|25[)). between two dissimilar net-neutral dielectric slabs interacting across 
a medium of dielectric constant e m varying in the range e m = 10,15,25,40 (dashed curves from bottom). Here we have 
fixed the dielectric constant of the slabs as ei = 5, £2 = 50, which contain uncorrelated quenched disorder (£1 =£2=0) of 
fixed disorder variances g\b = 32b = 5 x 10 -8 nm~ 3 and g\ s = <?2s = 0. Inset shows a closer view of the region around the 
minimum, b) Same as (a) but plotted for correlated quenched disorder. Here we fix e m = 20 and vary the correlation length as 
£/Ib = 0, 20, 100, 500, 10 4 , which is taken to be equal in both slabs £1 = £2 = £• Inset again shows a closer view of the region 
around the minimum, c) Same as (a) but for (uncorrelated) annealed disorder obtained from Eq. (|30[) for e m = 10, 20, 40. Inset 
shows the ratio of the total force to the pure zero-frequency vdW force (|18p in the absence of charge disorder for a wider range 
of separations. The horizontal dotted lines show the limiting expression (|3T[) . 



possibility of repulsive total forces in the annealed case 
that will be discussed in the following section. 



Non-monotonic interaction between dissimilar 
slabs with ei < e m < £2 



It is clear that the vdW force (|Po)) becomes repulsive 
when the dielectric constant of the intervening medium 
is between that of the two semi-infinite slabs, i.e., £1 < 
£m < £2 (such a case has been investigated experimen- 
tally in Ref. |l5|). However, in this case, the force stem- 
ming from the quenched charge disorder can be attractive 
and thus, if present, may compete against the repulsion 
due to the vdW forces. In order to investigate this situ- 
ation further, let us first consider the case of an uncorre- 
lated disorder by setting the disorder correlation lengths 
equal to zero. Also without loss of generality, we assume 
that the slabs contain only quenched bulk disorder with 
equal variances gu, = 92b = 9b in the two slabs. The net 
quenched force in this case reads 



/3/quonchod _ ^h(^1^2) 



9b^BX 



S 8nD 3 2(ei+e 2 ) J D' 

where we have defined 

(ei — e m ) (£2 — £m) 



X 



62 



61 



(43) 



(44) 



Obviously, the contribution from bulk disorder (second 
term in Eq. (l4"3"l) ) can change sign as the dielectric con- 
stant 6 m is varied in the range 6% < e m < 62, while the 
vdW term (first term) remains always repulsive. It is 



easy to see that for e m — > 6\, we have \ > and thus 
an attractive disorder force, and for e m — > £2, we have 
X < and thus a repulsive disorder force. This behavior 
is shown in Fig. [7^,, where £1 = 5, £ 2 = 50 are fixed 
and £ m varies in the range e m = 10,15,25,40 (dashed 
curves from bottom). In accordance with our findings in 
the symmetric case, the large distance behavior is always 
dominated by the disorder contribution (as the disorder 
force decays more weakly with the separation) , while the 
repulsive vdW force in this case plays the role of a stabi- 
lizing force at small separations. 

The resulting effect is that the total force varies non- 
monotonically and vanishes at a finite distance, Dq, be- 
tween the two slabs given in the present case by 



n 2 
^0 



Li 3 (AiA 2 ; 
^9b^B X 



(El +£ 2 )- 



(45) 



This represents a stable 'equilibrium' separation (bound 
state) between the two slabs corresponding to a mini- 
mum in the interaction free energy. On the other hand, 
the maximum attractive force due to the influence of 
quenched disorder is reached at a larger separation, 



VSD . 



(46) 



These results may be used to optimize the thickness of 
the intervening medium in order to achieve the maximum 
or minimum force magnitude between the slabs. 

In the presence of a correlated disorder, the disorder- 
induced effects weaken (Section IIV A[) and thus, as the 
disorder correlation is increased, the attractive tail as 
well as the stable bound state are gradually washed out 
as shown in Fig. [TJa. 
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FIG. 8: a) The rescaled total force, f3fl%/S (Eq. (J25J) ) , between two dissimilar net-neutral dielectric slabs interacting across a 
medium of dielectric constant e m varying in the range e m = 30,40,60, 100 (dashed curves from bottom). Here we have fixed 
the dielectric constant of the slabs as ei = 15, £2 = 25, which contain uncorrelated quenched disorder (£1 = £2 = 0) of fixed 
disorder variances gib = 92b = 5 x 10 -8 nm -3 and gis = 52s = 0. Inset shows a closer view of the region around the maximum, 
b) Same as (a) but for annealed disorder obtained from Eq. (|30[) . Inset shows the ratio of the total force to the vdW force 
(|18[) for a wider range of separations, c) Same as (a) but for one slab being disorder free (gib = gis = 0, £1 = 0) and the other 
slab containing annealed charge disorder of variances g2b — 5 x 10 -8 nm -3 and gis = (dashed curves from top correspond to 
e m = 30, 40, 60, 100). Inset shows the ratio of the total force to the vdW force (|18[) for a wider range of separations. 



We find similar features when the disorder charges are 
annealed as seen in Fig. [7k. However, in line with our 
finding in the symmetric case, the net force falls off more 
rapidly with the separation when the disorder charges 
are annealed. The large distance behavior coincides again 
with the universal expression (|37l) and can be much larger 
in magnitude than the pure vdW force (|18[) when e m 
tends to the larger dielectric constant £2 as shown in the 
inset of Fig. [7b. 



C. Non-monotonic interaction between dissimilar 
slabs with ei,£2 < £m 

In the case where the intervening medium has a higher 
dielectric constant than the two slabs, £1,62 < £ ra , the 
vdW effect leads to an attractive force between the two 
slabs, Eq. ([18]) , which again dominates at small sepa- 
rations. If the disorder is quenched, the force generated 
by the disorder can become repulsive and thus lead to a 
potential barrier and a long repulsive tail at large sep- 
arations. This is shown in Fig. [8^, for the case of two 
slabs with uncorrelated bulk disorder of equal variances 
(dib = 92b = 9b), where we have fixed £\ — 15, £2 = 25 
and vary e m in the range e m = 30, 40, 60, 100 (dashed 
curves from bottom). The total force in this case is 
given again by Eq. (|43|) . Note that the disorder-induced 
force (second term) is always repulsive in the case with 
El, £2 < £ m (i-e., X < 0). Thus, the separation distance 
Do, Eq. (|45"|) . corresponds to an unstable 'equilibrium' 
distance between the two slabs and D max gives the dis- 
tance at which the maximum repulsive force due to the 
influence of quenched disorder is achieved. 



The above features change dramatically if the disor- 
der charges are assumed to be annealed. In fact, the net 
annealed force appears to follow a trend similar to what 
one expects from the pure vdW force, i.e., in contrast 
with the quenched case, the net annealed force turns out 
to be attractive and vary monotonically with the sepa- 
ration (Fig. [8)3). However, the ratio of the net force to 
the vdW force (fT8| shows that the relative magnitude of 
the force can vary non-monotonically and deviate signifi- 
cantly from the underlying vdW force (Fig. [5}d, inset). In 
particular, one observes that the relative net force may 
be enhanced here by an order of magnitude at large sep- 
arations. This is again due to fluctuations of annealed 
charges that can redistribute in the slabs in such a way 
as to minimize the free energy and thus favor a higher at- 
tractive force. This effect appears to be much stronger in 
an asymmetric system than in a symmetric system (Figs. 
[3]and[H]), where the net force remains of the same order 
as the pure vdW force. 



In the preceding discussion, we assumed that the two 
slabs have similar annealed disorder variances, which 
leads to a purely attractive force between the slabs for 
£i,£2 < s m . It turns out that the annealed disorder can 
also lead to a repulsive force in this case provided that 
the charge is distributed asymmetrically between the two 
slabs. In Fig. [5J;, we show the results for the case where 
one of the slabs is disorder free, but the other slab con- 
tains annealed charge disorder of finite variance. As seen, 
one can achieve an interaction (potential) barrier even 
with the annealed charges. This effect is stronger when 
£ m is smallest and disappears when e m — > 00. 
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VI. DISCUSSION 

In this paper, we have shown that the effective inter- 
action induced by the quenched or annealed monopo- 
lar charge disorder can give rise to various novel fea- 
tures in the overall interaction between two net-neutral 
semi-infinite dielectric slabs depending on the detailed 
assumptions about the disorder and the dielectric inho- 
mogeneities in the system. The quenched case and the 
annealed case of disorder differ in the sense that in the 
former the disordered charges are frozen and can not fluc- 
tuate, while in the latter the disordered charges are sub- 
ject to thermal fluctuations and adapt themselves to min- 
imize the free energy of the system. Our analysis is based 
on recent developments which unify the salient features 
of the zero-frequency vdW interaction and the physics of 
the charge disorder @, [l2|, . Based on these develop- 
ments we argue that, in the case of two dissimilar slabs, 
the charge disorder-induced electrostatic interaction can 
have an opposite sign to the zero-frequecny vdW interac- 
tion and can thus give rise to a non-monotonic net inter- 
action between the slabs. The most surprising features 
of our analysis pertain to the case of disorder-induced in- 
teractions across a medium of higher dielectric constant 
than that of the two slabs. In this case the net force may 
become strongly repulsive and can lead to a potential 
barrier for stiction when combined with the attractive 
vdW force. In the opposite case, where the intervening 
medium has a dielectric constant in between that of the 
two slabs, the disorder may generate a long-range attrac- 
tion, which opposes the repulsive vdW force and can thus 
promote a stable bound state for the two bounding di- 
electrics. These salient features of the disorder-induced 
interactions are due to the slower decay of these inter- 
actions and their stronger dependence on the dielectric 
inhomogeneities in the system as compared to the corre- 
sponding vdW interaction. 

This comparison between the disorder-induced and the 
classical zero-frequency vdW interaction is in order since 
both of them are expected to be valid in the regime of 
large intersurface separations (or high temperatures) as 
considered in this work and which is also most relevant to 
recent experiments 0. The precise correction presented 
by the higher-order Matsubara frequencies [2l[ is highly 
material specific, but its magnitude (relative to the zero- 
frequency term) is typically small for the most part of the 
separation range considered here and remains negligible 
in comparison with the disorder effects. 

Being disorder induced one should be in principle also 



able to compute the corresponding statistical moments 
of the disorder interactions averaged over the disorder. 
Note that these fluctuations over the disorder will be dif- 
ferent from the thermal fluctuations of the Casimir force 
studied in Ref. [22] which are present even in the absence 
of disorder. 
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Appendix A: General expression for the force in the 
quenched case 

In the general case where the two slabs have different 
dielectric constants and charge disorder parameters, we 
can write the total force, Eq. (|25l) , as a series expansion 
in powers of £i a /D (where for the two slabs, i = 1,2, and 
bulk and surface disorder, a = b, s), i.e. 



,/quenchcd 

/( 0) + £ (/ <») + 



in) , ,(n) ); (A1) 



n=l 



where we obtain 

516^b( £ 2 + £m)A 2 gis^BSm l°g(|l - AlA 2 | 



S 2(e 1 +e m )(ei+e 2 ) J D ( El + e m ) 2 Ai£> 2 

.92fcMei + £ m)Ai g 2 ^se m log(|l - AiA 2 |) 



2(e 2 +e m )(ei + E 2 )D (e a + e m ) 2 A 2 L> 2 

and 



(A2) 



00 



E(- 1 )"-^r C2 »+i Li2 »( A i A2 ) 



S (e 1 +e m ) 2 A 1J D^ v 

v ; n—1 



D 



4:gi s £ B £r 



E(- 1 ) n & C2 «+ 2Li2 »+i(AiA 2 ). (A3) 



(e 1 +e m ) 2 A 1 ^ 2 ^ ' D* 

The expression for is obtained simply by replacing 
the subindex 1 with 2 and vice versa. 
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